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Abstract 

We present a study of the noise driven escape of an overdamped Brownian particle 
moving in a cubic potential profile with a metastable state. We analyze the role of 
the initial conditions of the particle on the enhancement of the average escape time 
as a function of the noise intensity for fixed and fluctuating potentials. We observe 
the noise enhanced stability effect for all the initial unstable states investigated. For 
a fixed potential we find a peculiar initial condition x c which separates the set of 
the initial unstable states in two regions: those which give rise to divergences from 
those which show nonmonotonic behavior of the average escape time. For fluctuating 
potential at this particular initial condition and for low noise intensity we find large 
fluctuations of the average escape time. 
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1 Introduction 



In the two last decades a remarkable increase has been done in the study of 
noise induced effects on nonlinear nonequilibrium systems, with attention also 
to complex and biological systems [1,2,3,4]. The introduction of noise and of 
a deterministic driving force in such a systems gives rise to counterintuitive 
effects which explain apparently anomalous behaviors in experiments. Nonlin- 
ear relaxation decay of physical systems from an initial unstable or metastable 
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state involves fundamental aspects of non-equilibrium statistical mechanics. 
Examples of these resonance-like phenomena are stochastic resonance (SR) 
[3,5], resonant activation (RA) [6], noise enhanced stability (NES) [7,8], noise 
induced phase transitions (see ref.[l]), etc... A nonmonotonic behavior of the 
average escape time as a function of the noise intensity was revealed in a nu- 
merical study of a Brownian particle moving in a periodic fluctuating cubic 
potential [9]. This is the NES phenomenon: the stability of an otherwise un- 
stable system can be enhanced by the presence of a finite amount of noise. 
The nonmonotonic behavior that often occurs in NES effect contradicts stan- 
dard Kramers-like behavior, i.e. an exponential or monotonic decrease of the 
mean escape time with noise intensity [10]. The NES phenomenon was exper- 
imentally detected in the transient dynamics of an unstable physical system 
[7] and observed in different physical systems [11]. Recently some works have 
been done studying piecewise linear potentials with a metastable state, and 
exact evaluation of the decay time was obtained [8]. In this work we present a 
study of the average decay time of an overdamped Brownian particle subject 
to a cubic potential with a metastable state at different unstable initial con- 
ditions. A relevant result of this work is the existence, for fixed potential, of a 
peculiar initial position x c of the Brownian particle, which corresponds to the 
intersection point between the potential profile and the x-axis (see Fig. 1). For 
all the initial unstable states between the maximum of the potential and x c , 
the average escape time diverges, while for initial positions x Q > x c we find a 
nonmonotonic behavior of the same quantity. Our results are consistent with 
that obtained in the case of piecewise linear potential profile [8]. For a periodic 
fluctuating potential we recover the NES phenomenon for different values of 
the amplitude and of the frequency of the periodical driving force and large 
fluctuations of the average escape time as a function of the amplitude and the 
frequency at x Q = x c . 



2 The Model 



The starting point of our study is the following Langevin equation 



where is the white Gaussian noise with the usual statistical properties: 
(£(*)} = and + r)> = 5(r), and 

U(x,t) = ax 2 — bx 3 — xAcos(uit), (2) 
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Fig. 1. The cubic potential U(x) with the various initial positions investigated (dots); 
x c is the critical initial position, (a) and (b) indicate the limit curves of the oscillating 
potential. The absorbing boundary is xf = 20. 

is the potential shown in Fig. 1, with A and u respectively the amplitude and 
the frequency of the driving force, and a = 0.3, b = 0.2. The potential profile 
has a local stable state at Xq = and an unstable state at Xq — 1. For fixed 
potential A = the average escape time of a particle starting from x Q and 
reaching a final position xf is given by [12] 

X F z 

t(x ,x f ) = pJ eMz) I e~ My) d y d z- (3) 
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where u(y) = {0.3y 2 — 0.2y 3 )/D is a dimensionless potential profile. We eval- 
uate this double integral partly analytically and partly numerically, splitting 
it as 



T[X 



xp z 

x F ) = ^f e 2u{z) dz J e~ 2u M + J 
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dy 



(4) 



and evaluating analytically the first term inside the square parenthesis: 

/ e~ 2u ^ = 0.6046 e- z [I. 1/3 (z) + I 1/3 (z)} - \ 2 F 2 (1 1; | i -2z), (5) 

— OO 



where z = 1/(10D), I n (z) is the modified Bessel function of the first kind and 
pF g (ai, a<i \ bi, b 2 ; z) is the generalized hypergeo metric function. 
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Fig. 2. a) Mean First Passage Time as a function of noise intensity for all the nine 
initial positions investigated (see Fig. 1), namely: x Q = 1.1 -f- 1.9, with steps of 0.1. 
The number of realizations is 2000 and xf = 20. b) The standard deviation of the 
first passage time distribution, which shows a nonmonotonic behavior similar to 
that of MFPT. 



3 Results and Comments 



For all the initial unstable states beyond the potential barrier (see Fig. 1) we 
find an enhancement of the average escape time with respect to the deter- 
ministic time as a function of the noise intensity (Fig. 2). Specifically, for high 
values of the noise intensity with respect to the height of the barrier, we recover 
the Kramers behavior. For intermediate values of noise intensity, we observe 
an increase of the MFPT whose maximum goes to infinity when the initial 
position is in the range x max < x Q < x c . These results are quite consistent with 
those obtained in the case of piecewise linear potential profiles [8], where a 
divergent behavior for D — > is found. On the other hand for very low values 
of noise intensity we obtain the deterministic decay time. This effect, which 
appears in contradiction with the divergence for D — > 0, can be explained 
with the finite ensemble of particles in numerical experiments. Indeed, for 
low noise intensity the probability that some particle is trapped is very low 
and decreases exponentially to zero, that is only very few particles will be 
trapped. Only these particles, which represent very rare events for very small 
noise values, contribute to the enhancement of the escape time and give rise 
to the divergence. Therefore, because of this, we do not observe such particles 
in simulation and the average escape time becomes equal to the deterministic 
time. In Fig. 2 we can clearly see this effect. For initial positions approach- 
ing the maximum of the potential this effect is less pronounced, because the 
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trapping probability of the particle is proportional to exp(— Au(x )/D), where 
Au(x Q ) is the potential barrier "seen" from the particle at x(0) = x Q . When 
the particle is near the maximum a relatively low amount of noise can push 
the particle back into the "stable" state, where it remains trapped for a long 
time, because of the low noise intensity. For higher values of initial position 
the amount of noise to put the particle into the potential well must be higher. 
This explains why the maximum of the average decay time is shifted towards 
higher values of the noise intensity. The observation time in digital simulation 
is finite, therefore the simulated escape time takes the maximum at this limit- 
ing point, and this is what we find (see Fig. 2a). For all the initial conditions 
within the range: x m < x Q < x c , we find that the average escape time can be 
more and more increased by the noise, while for initial conditions x D > x c we 
obtain nonmonotonic behavior with a maximum and a finite average escape 
time when the noise intensity goes to zero. We calculate also the standard 
deviation (cr T ) of MFPT as a function of the noise intensity for all the initial 
positions investigated (see Fig.2b). For x Q > x c we obtain a nonmonotonic 
behavior like the MFPT behavior. For x Q < x c the standard deviation grows 
towards divergent values for D — > oo, and this corresponds to a big tail in 
the first passage time distribution. From the above analysis it appears that 
a peculiar initial position for the potential profile of Eq.(2) is x Q = x c = 1.5. 
We calculate then the theoretical average escape time from Eqs.(4,5) for an 
initial position just after the cross point x c and we show the result in Fig. 3. 
The agreement between the numerical simulation of the SDE of the Eq. (1) 
and the theoretical evaluation is very good. 
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Fig. 3. Average escape time as a function of the noise intensity D, for an initial 
position x Q = 1.51 > x c = 1.5 and xp = 2.2. Comparison between theoretical 
MFPT (full circles) and the simulation of the SDE (empty circles)(Eq. (1)). 

We calculate the average escape time at the static cross position x c = 1.5 for 
a periodical fluctuating cubic potential and the results are shown in Fig.4 and 
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Fig. 5. Here we recover the enhancement of stability of the metastable state for 
all values of the amplitude and of the frequency of the periodical driving force 
investigated, with a nonmonotonic behavior of MFPT (a- figures). We find also 
that (b-figures): (i) the MFPT values are almost independent from the ampli- 
tudes and the frequencies of the driving force, except for the case D = 0.05; 
(ii) very large fluctuations of MFPT as a function of both the amplitude and 
the frequency of the periodic driving force at D = 0.05 are present. The ori- 
gin of these large fluctuations can be ascribed to the particular initial position 
represented by x c , which is one of the two boundaries delimiting the parameter 
region (A,u) where NES effect can be observed [8]. Because of the oscillation 
of the potential between two limiting curves the intersection point is a func- 
tion of the time x c (t) and as a consequence the Brownian particle experiences, 
during each period of the driving force, the two different dynamical regimes 
of the case of the static potential discussed above for D — > 0. Specifically the 
escape time passes from a divergent regime (x c (t) < x c ) to a convergent one 
(x c (t) > x c ), where x c is the intersection point for the static potential. This 
effect becomes more important for low noise intensity, i.e. for noise intensities 
less than the barrier height, as we can see in Figs. 4 and 5 for D = 0.05. We 
obtain the same large fluctuations for other values of D less than the barrier 
height (0.1), not reported here for sake of clearness of the figures. For values 
of amplitude and frequencies of the driving force higher then those presented 
in the figures, up to v = 20 and up to A = 3 respectively, the calculations 
confirm these large fluctuations. We also investigated the influence of different 
initial phases on the MFPT behavior as a function of the frequency and we 
recover large fluctuations again for D = 0.05. 



4 Conclusions 



Nonmonotonic behavior of the mean escape time as a function of noise in- 
tensity is a noise-induced effect for nonlinear nonequilibrium systems with 
metastable states. In this work we analyzed the role of the initial conditions 
on the enhancement of the escape time from initial unstable states for a cubic 
potential. We obtain NES effect for static and periodical fluctuating poten- 
tial and an enhancement of the NES effect for initial positions between the 
maximum of the potential well and the cross point x c . We find also large 
fluctuations of the average escape time as a function of both the amplitude 
and the frequency of the sinusoidal force for low noise intensity, due to the 
different dynamical regimes experienced by the Brownian particle. Our results 
obtained for a particle moving in a cubic potential are quite general, because 
we always obtain NES effect when a particle is initially located just to the 
right of a local potential maximum, with a local minimum in its left side and 
the global escape region in its right side. 



6 




Fig. 4. MFPT as a function of noise 
intensity (a) and of the amplitude (b) 
of the driving force. The parameter 
values are: x Q = 1.5, v = 10Hz. 
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Fig. 5. MFPT as a function of noise 
intensity (a) and of the frequency (b) 
of the driving force. The parameter 
values are: x Q = 1.5, A = 0.02. 
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